Ts = [12,15,20,23,26,30,35,40,50,60,80,100,130,160,190,220,250,275,296];
figure
hold on
dataall=zeros(19,1997);
for num=1:length(Ts)
    filename=['HV',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1)/8.06554,data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1)/8.06554;
    dataall(num,:)=data(:,2);
    ylim([20,500])
    title('HV')
end
figure
hold on
dataall=zeros(19,1997);
for num=1:length(Ts)
    filename=['VV',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1)/8.06554,data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1)/8.06554;
    dataall(num,:)=data(:,2);
    ylim([20,500])
    title('VV')
end
figure
imagesc(Es,Ts,dataall)
colormap('jet')
clim([20,40])

figure


[X,Y]=meshgrid(Es,Ts)
mesh(X,Y,dataall)
colormap('jet')
clim([15,60])
xlim([55,70])